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ABSTRACT 

We consider the possibility of using the hierarchical approximation to under- 
stand the continuum limit of a reformulation of the 3D Ising model initiated by 
Polyakov. We introduce several new formulations of the hierarchical model us- 
ing dual or fermionic variables. We discuss several aspects of the renormalization 
group transformation in terms of these new variables. We mention a reformulation 
of the model closely related to string models proposed by Zabrodin. 



1. Introduction 



There has been a long-standing effort to describe the large distance behavior of 
gauge theories in terms of a string theory.' 11 This is obviously not an easy task and it 
would be very interesting to have at hand a simple example where all details can be 
carried out explicitly. It seems worth trying with a 3D Z(2) gauge theory, because 
this model is dual to the nearest neighbor 3D Ising model for which we have a 
reasonable understanding of the critical behavior. Polyakov has shown that there 
exist linear relations among the averages of products of order and disorder variables 
taken along contours. These reveal the existence of fermionic excitations "held 
together" by gauge-invariance. Subsequently, these contour averages have been 
reexpressed as sums over "equipped surfaces" by Polyakov and Dotsenko.' 41 These 
results strongly suggest that near the critical temperature, the 3D Ising model can 
be described in terms of some free fermionic string. 

Despite significant efforts, we do not yet have a clear understanding of the con- 
tinuum limit of the reformulation of Polyakov and Dotsenko. For this reason, we 
have suggested 161 considering this problem in the context of the hierarchical approx- 
imation. This approximation has been justified by Wilson 7 and used to calculate 
the critical exponents of the 3D Ising model with good precision. The practi- 
cal advantage of the hierarchical approximation is that the renormalization group 
procedure reduces to the study of a simple integral equation ("the approximate 
recursion formula"). The fixed points of this recursion formula and the relevant 
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eigenvectors of the linearized transformation have been studied in great detail. The 
critical exponents have a very simple form in terms of the corresponding eigenval- 
ues. 

The hierarchical approximation holds exactly in the case of the hierarchical 
model 191 (HM in the following) which will be the main object of the present study. 
For definiteness, this model is briefly reviewed in section 2. We then introduce 
several new formulations of the HM. The model dual to the HM is described in 
section 3. Various fermionizations are discussed in section 4. Due to the fact that 
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the model has long range interactions, specific methods were necessary to obtain 
these results. These methods indeed apply to a larger class of models and are 
reported with due mathematical care in a separate article.' 101 Finally, we mention 
in the conclusion that that the model can be rewritten as a nearest neighbors 
model on a "branch" by introducing auxiliary (gaussian) fields. This construction 
is related, but distinct on several points, to string models proposed by Zabrodin. [nl 

The main question we address is how the extraordinary simplification provided 
by the renormalization group method (we only need the relevant directions) takes 
place when, instead of the order variables, we use the new variables discussed above. 
In all the cases, the renormalization group transformation appears to be more 
involved than in the usual formulation. This can be understood from the fact that 
in the HM, the number of links grows proportionally to the square of the number 
of sites. However, preliminary results indicate that when the symmetries of the 
model are taken into account, the procedure can be reduced to a size comparable 
to the usual one. The hard work (fixed points, linearization) remains to be done. 
A proper understanding of this question may shed a new light on the continuum 
limit of the sum over equipped surfaces discussed above. 



For definiteness recall a few facts concerning the HM. More details and refer- 
ences can be found in the lecture notes of Collet and Eckmann or Gawedzki and 
Kupiainen. The hamiltonian of a hierarchical model with q n sites can be written 

as 



All the indices ij run from 1 to q, an integer which controls the size of the boxes used 
during the renormalization group transformation. We assume that 1 < c < q. It is 
suggestive to write c = q l ~ dh where is the Hausdorff dimension' 121 of the random 
walk associated with H. If the spins are integrated with a gaussian measure, 



2. The Hierarchical Model 




(1) 



3 



the model has the same critical exponents as an ordinary gaussian model in D 
dimensions, provided that dh — A non-gaussian continuum limit exists for 
1 > d h > 1/2. 

In the next two sections, we use the short notation % instead of (i n , 
and we define a function v(i,j) which is equal to I if i\ is the first index (starting 
from the left) differing in % and j. Another way of saying it that the smallest box 
containing % and j contains q v ( l >fi sites. The partition function reads 



Wi=±i} 



with 



Ki i =P(i-^r 1 (& iJ) -(^) n+1 ). (3) 



q 2 q 2 q 2 ' 



The notation i < j in the sum means that i is distinct from j and that the pairs 
of distinct elements which can be obtained by interchange are only counted once. 



3. The Duality Transformation 

As we have seen, the Kramers- Wannier duality^ 1 plays an important role in 
Polyakov's approach. In the case of the HM, a dual spin S^jja is associated with 
any triangle ijk. These triangles can be used to build surfaces whose boundary 
is a contour appearing in the high temperature expansion of the model. The fact 
that tetrahedrons are closed surfaces provides the principle underlying the gauge- 
invariance of the dual formulation. All we need to check is that all the terms of 
the high temperature expansion appear with the same "gauge-multiplicity" . This 
is proven elsewhere! 101 In the following, we only state the main results. 

If we define the dual couplings tanhKij = e~ 2Dij , we can write the partition 
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function of Eq.(2) as 



Z = ( JJ coshKij) 2 N ~( 3 J ^ e fc ^ . (4) 

»<j {s !3fc =±i} 

This reformulation has a local invariance. For any four distinct sites k and I, 
the dual interaction energy is invariant under the simultaneous changes of sign of 
S(ijk), s (ijl)i s (ikl) and s (jkl)- A possible gauge-fixing condition is = if none 

of i,j and k are equal to a given site. In general, if we start with q n sites, the dual 
formulation has (1/2) (q n — l)(q n — 2) physical degrees of freedom. This number 
can be obtained directly by subtracting the q n — 1 independent constraints on the 
parity of the number of visits at each site from the total number of links. 

We insist on the fact that the symmetries of the model [91 are crucial to deal with 
this proliferation. In order to fix the ideas, let us just compute the high temper- 
ature expansion of the partition function "hierarchically". We perform successive 
integrations over the links variables inside boxes of size q l , starting with / = 1 
and proceeding similarly for higher I. At first sight it looks like we have to keep 
track of a prohibitively large number of quantities, namely 2 q , since each site can 
be visited an even or an odd number of times and we need this information when 
we integrate over the links in larger boxes. At the end of the calculation (say for 
I = n), we only retain those contributions where the number of visits at each site is 
even. This indeed amounts to calculate all the correlation functions for q l sites for 
I = l,2..,n — 1. Fortunately, when the symmetries of the model are taken into ac- 
count, the number of correlation functions to calculate is reduced logarithmically. 
It seems thus, in principle possible to use the renormalization group procedure in 
the dual formulation. Nevertheless, for practical reasons, it seems advantageous 
to first represent the high temperature expansion as an integral over Grassmann 
variables 1131 . 
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4. Hierarchical Fermions 



A compact formulation of the high temperature expansion of the HM can be ob- 
tained [u>1 

by associating to each link ij the four Grassmann variables ip^ippXi an d Xj- 
The final result for the partition function is 

Z=C[[ coshK^ J [d^dx] e J] *( E (5) 

where [ci^cix] — EI d-xldipjdXjdipj and is equal to sinh(x) (resp. cosh(x)) if 

g is even (resp. odd). The average value of products of order or disorder variables 
and the Schwinger-Dyson equations are easily obtainable from Eq.(5). Note that, 
unlike in the 3D Ising model, these equations are highly non-linear. Note also that 
the number of Grassmann variables introduced is not minimal. We have doubled 
the number of fermions in order to avoid the splitting of the partition function into 
an uncontrollably large number of terms. The details concerning the above results 
can be found in Ref.[10]. 

This representation seems well-suited for the renormalization group approach. 
We can proceed in two steps. First, we integrate over the variables corresponding 
to the links inside the boxes of size q. Second, we make a linear transformation 
among the 4g 2 variables associated to the links joining any pair of boxes of size q 
in such a way that the average over upper and lower indices for the four type of 
variables are among the new variables. We then perform a gaussian integration 
over the remaining ones. After an appropriate rescaling, the "kinetic term" keeps 
its original form. What is transformed is the function z(x) appearing in Eq.(5). 
We intend to study this transformation with the methods used for the recursion 
formula mentioned in the introduction. 

Note also that we can construct a fermionic representation of the low temper- 
ature expansion. To each dual site ijk we associate the three Grassmann variables 



6 



^iji ^iifc' ^jk anc ^ their \ counterpart. The kinetic term contains a part of the form 

Y^thD^ n ^ (6) 

i<j k;k^i,j 

and the method based on gaussian integration is not applicable. 

5. Conclusions and Perspectives 

We have introduced several reformulations of the hierarchical model. The 
renormalization group method seems an appropriate tool to handle at least one 
of them. We think that a study of the fixed points and the relevant directions of 
the new transformation will provide a better understanding of the continuum limit 
of the reformulation of the 3D Ising model proposed by Polyakov and Dotsenko. 
This issue may also be clarified by a more systematic understanding of the hierar- 
chical approximation. This could be done by using a complete set of multiplicative 
characters and taking the degree of ramification as the order of perturbation. 

Note also that the HM can be reformulated as a nearest-neighbor model on a 
"branch". The hamiltonian reads 

n i li 

n - 2^ ^ 2 [a (in,....,ii+i) ) q a (in,-..,ii+i) 2^ (T (i n ,....,ii) ) l j 

1=1 i n ,...,ii +1 it 

with a\ — 1 and a\ — 1+| for I > 2. The af^ are the spins variables of the HM as in 

section 2. The new variables a^ l \ for / > 2 are integrated with (27T/9) -1 / 2 f da\ . 

This construction has been inspired by Zabrodin's models' 111 designed to reproduce 
the p-adic generalization 1141 of the tachyon amplitudes of string theory. Our refor- 
mulation differs on several points (discreteness of the boundary, the local couplings 
are not necessarily the same) from these models, however, it might be worth ex- 
ploiting their resemblances. The "vertex operator" Vk = ^2 Exp(kaf^) can be 
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used to generate the average value of functions of the total spin. Note however, 
that in the gaussian case, the SL2 invariance 114111 applies only for d^—1 i.e D = 2. 
Otherwise, a "mass term" spoils the transformation properties and we cannot use 
the Koba-Nielsen trick. 
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